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REMARKS 

Applicants acknowledge the indication of the allowability of the subject 
matter of claim 7, as set forth in Item 7 on page 6 of the Office Action. In 
particular, claim 7 would be allowable if rewritten in independent form. 
Nevertheless, for the reasons set forth hereinafter, applicant respectfully 
submits that claim 7 is allowable in its present dependent form. 

Applicant also takes note of the Examiner's request that it provide a copy 
of the foreign patent documents cited in the Information Disclosure Statement 
filed August 19, 2004. Of these, one is the Snell reference (GB 2 268 138A), 
which is cited and relied on in the Office Action. A copy of International Patent 
Document WO 98/49058 is enclosed for the Examiner's convenience. Applicant 
notes in this regard, however, that, as pointed out in the August 19, 2004 IDS, a 
copy of this reference was provided to the International Searching Authority and 
is therefore not required to be submitted with the IDS itself. Therefore, 
applicant respectfully submits that it is entitled to have International Patent 
Document WO 98/49058 considered on the record herein, and respectfully 
requests that the Examiner provide an additional copy of the PTO Form 1449, 
initialed to show that this has been done. 

In response to the rejection of claims 9 and 11, as set forth in Item 5 on 
pages 2 and 3 of the Office Action, applicant has canceled claims 9 and 11. 
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Accordingly, this ground of rejection, as well as the rejection of the latter claims 
based on the Carrico et al. reference (U.S. Patent No. 6,937,968), are rendered 
moot. 

Claims 1, 2, 5 and 10 have been rejected under 35 USC § 102(e) as 
anticipated by Snell (GB 2 268 138 A), while claim 3 has been rejected under 35 
USC § 103(a) as unpatentable over Snell in view of Madhani et al. (U.S. Patent 
No. 6,786,896) and claims 4 and 6 have been rejected as unpatentable over Snell 
in view of Ocampo (U.S. Patent No. 6,341,749). However, for the reasons set 
forth hereinafter, applicants respectfully submit that all claims which remain of 
record in this application, including new claims 13 and 14, distinguish over the 
cited references, whether considered separately or in combination. 

By the foregoing amendment, claim 3 has been canceled, and its 
limitations incorporated into independent claims 1 and 10, as well as new 
independent claims 13 and 14. With regard to claim 3, Item 5 of the Office 
Action indicates that it would have been obvious to a person skilled in the art to 
use null space vectors of Madhani et al. in mathematical computations for 
control balancing of a spacecraft as taught in Snell. In regard to this premise, 
applicant notes, however, that the Snell reference is directed to thrust vectoring 
of a multi-engined VTOL aircraft (and, specifically, to the problem of a possible 
failure of an engine of a VTOL aircraft while it is operating in a vertical mode), 
while Madhani et al. relates to a method for performing surgery, and to a 
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servomechanical surgical apparatus for holding and manipulating human tissue 
under control of a teleoperator system. Applicants respectfully submit that it 
would not be obvious to a person skilled in the field of thrust vectoring for the 
control of spacecraft to look to the field of control of robotic medical apparatus to 
find a solution to the problem addressed by the present application — that is, the 
difficulty that the constraints applicable to the prior art cannot be satisfied for 
high control demands with respect to the ideal maximal thrust which can be 
generated from the thrusters in a spacecraft. Moreover, even were such a skilled 
person to look to the field of robotic controls for this purpose, applicant 
respectfully submits that there is no suggestion contained in either document 
(Snell or Madhani et al.) that they should or could be combined, or if possible 
how such a combination would be effected. In the absence of any such teaching, 
applicant respectfully submits that the combination of Snell and Madhani et al. 
does not render the invention of independent claims 1, 10, 13 and 14 obvious, 
and that the latter claims are therefore allowable. 

If there are any questions regarding this amendment or the application in 
general, a telephone call to the undersigned would be appreciated since this 
should expedite the prosecution of the application for all concerned. 
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If necessary to effect a timely response, this paper should be considered as 
a petition for an Extension of Time sufficient to effect a timely response, and 
please charge any deficiency in fees or credit any overpayments to Deposit 



Account No. 05-1323 (Docket #101806.55264US). 



CROWELL & MORING LLP 
Intellectual Property Group 
P.O. Box 14300 
Washington, DC 20044-4300 
Telephone No.: (202) 624-2500 
Facsimile No.: (202) 628-8844 
GRE:aw 

2754876 



Respectfully submitted, 



May 19, 2006 




Gary R. Edwards 
Registration No. 31,824 
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(54) Title: METHOD FOR LOW-FUEL, COMPUTER-ASSISTED CONTROL OF ANY NUMBER OF THRUSTERS ARRANGED ON 
A SPACECRAFT 

(54) Bezeichnung: VERFAHREN ZUR TREIBSTOFFMINIMALEN, RECHNERGESTOTZTEN ANSTEUERUNG FOR BELIEBIG AN 
EINEM RAUMFAHRZEUG ANGEORDNETE DtlSEN 

(57) Abstract 

The invention relates to a method for low-fuel, computer-assisted control of any number of thrusts mounted in any given arrangement 
on a spacecraft, which solves the problem of linear optimization by reducing computational effort in a very short time. The method is based 
on the dual simplex algorithm and suitable provision of a dual admission start table. The dual admission start table is set up on the basis of 
a previously calculated optimum table and a transformation of the actual force/power vector pertaining to the chosen optimum table. Given 
processing steps of the dual simplex algorithm are skipped totally or are substituted by processing steps which require substantially less 
computational effort. Finally, a reduced table as opposed to the simplex table is used. Said reduced table consists solely of an identification 
number for selecting a thruster from the enure number of thrusters and the vector of the control values of the thrusters which are contained 
in said selection. Reduction in computational effort is achieved by carrying out given processing steps in advance and by providing the 
corresponding results in the form of data fields and index or identification number matrixes in the working memory of the computer. Said 
results from the initialization phase of the method need only to be addressed and retrieved in order to calculate thruster control. 

(57) Zusammenfassung 

Es soil ein Verfahren zur insgesamt treibstoffminimalen, rechnergestutzten Ansteuerung einer beliebig festgelegten Anzahl und beliebig 
festgelegten Anordnung von Dusen an einem Raumfahrzeug angegeben werden, das dieses lineare Optimierungsproblem durch Reduzierung 
des Rechenaufwands in mdglichst kurzer Zeit lost. Das Verfahren basiert auf dem dualen Simplexalgorithmus und der geeigneten 
Bereitstellung eines dual zulassigen Starttableaus. Das dual zulassige Starttableau wird aus einem vorab berechneten Optimaltableau 
und der zu diesem gewahlten Starttableau dazugehorigen Transformation des aktuellen Krafte-Momenten-Vektors aufgebaut. Bestimmte 
Verarbeitungssschritte des dualen Simplexalgorithmus entfallen vollstandig bzw. werden durch Verarbeitungsschritte substituiert, die 
wesentlich weniger Rechenaufwand bendtigen. Schliefilich wird nicht mit dem Simplextableau, sondem mit einem reduzierten Tableau 
gearbeitet. Dieses reduzierte Tableau besteht lediglich aus einer Identifikationsnummer fur eine Dttsenauswahl aus der Gesamtzahl der Dusen 
und dem Vektor der Ansteuerwerte fur diejenigen Dusen, die in dieser Dusenauswahl enthalten sind. Die Reduzierung des Rechenaufwands 
wird dadurch erreicht, daft bestimmte Verarbeitungsschritte vorab durchgefuhrt werden und die entsprechenden Ergebnisse dann in Form 
von Datenfeldem und Indizes- bzw. Identifikationsnummernmatrizen im Arbeitsspeicher des Rechners bereitliegen. Diese Ergebnisse des 
Initialisierungsteils des Verfahrens brauchen dann zur Berechnung der Dusenansteuening nur noch adressiert und abgerufen zu werden. 
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1 0 Verfahren zur treibstof fminimalen, rechnerqes tiitz ten Ansteue- 

rung fur beliebig an einem Raumfahrzeuq anqeordnete Diisen 



15 Die Erfindung bezieht sich auf ein Verfahren zur insgesamt 

* 

treibstof fminimalen , rechnerges tut z ten Ansteuerung einer 
beliebig festgelegten Anzahl und beliebig festgelegten 
Anordnung von Diisen an einem Raumf ahrzeug. Die Diisen konnen 
dabei diskret oder auch kontinuierlich ansteuerbar sein. 

20 

Fur die Ausubung von auBeren Kraften und Momenten an einem 
Raumfahrzeug ist der Einsatz von Diisen in vielen Situationen 
notwendig. Da der Einsatz der Diisen Treibstoff verbraucht und 
dieser unmittelbar EinfluB auf das Startgewicht und die 

25 Lebensdauer von Raumf ahrzeugen und insbesondere von Satelliten 

hat, kommt der Dusenans teuerung mit in der Surame minimalem 
Treibstof fverbrauch eine wichtige Bedeutung zu. Diisenanordnun- 
gen, die prinzipiell einen niedrigen Treibstof fverbrauch 
ermoglichen, eine hohere oder gleiche Redundanz bei kleinerer 

30 Diisenanzahl aufweisen und/oder durch die Nutzlast am Raumfahr- 

zeug beeinf lufit sind, sind oft unregelmaBige und schief winklige 
Diisenanordnungen. Aus diesen Griinden ist die Forderung nach der 
treibstof fminimalen Dusenans teuerung , insbesondere auch fiir 
beliebige ( schief winklige ) feste Diisenanordnungen an Raumf ahr- 

35 zeugen von besonderer Relevanz. 

Es ist bekannt, da£ sich das Problem der treibstof fminimalen 
Diisenansteuerung auf ein lineares Optimierungsproblem zuriick- 
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fiihren laBt, fur das der Simplexalgori thmus ein geeignetes und 
effizientes Losungsverf ahren darstellt (siehe zum Beispiel 
P„J. Wiktor "Minimum Control Authority Plot: A Tool for 
Designing Thruster Systems" In: Journal of Guidance, Control, 
and Dynamics, Vol. 17, No. 5, Sept. -Oct. 1994, pp. 4-5)- Der 
Einsatz des Simplexalgorithmus scheitert jedoch haufig an der 
benotigten Rechenzeit. So konnen oft die Echtzeitanf orderungen 
zum Beispiel fur den Einsatz des Simplexalgorithmus innerhalb 
des Regeltaktes des Lage- und Positionsregelkreises eines 
Raumf ahrzeuges nicht erfiillt werden, weil die besonders in 
hohen Orbits einsetzbaren Rechner den Simplexalgorithmus nicht 
ausreichend schnell abarbeiten konnen. 

Eine Vielzahl der sich heute im Einsatz befindenden Raumfahr- 
zeuge (Satelliten) besitzen eine Dusenanordnung, bei der jede 
Duse zum korperfesten Koordinatensystem des Raumf ahrzeugs zu 
einer Achse parallel und damit zu den zwei anderen Achsen 
orthogonal angeordnet ist. Diese orthogonale und 
achsenparallele Anordnung hat zur Folge, daB jede Duse nur 
Krafte in einem Richtungssinn einer Achsenrichtung und Momente 
nur in einem Drehsinn um die beiden anderen Achsen erzeugen , 
kann. Der Vorteil dieser Anordnung ist, daB die Steuerung eines 
jeden Preiheitsgrades des Raumf ahrzeuges als starrem Korper 
durch den Einsatz von jeweils nur zwei Dusen moglich ist, ohne 
andere Freihei tsgrade zu beeinf lussen • Die Ansteuerung der 
Diisen ist fur diese Dusenanordnung sehr einfach. Der Nachteil 
dieser Dusenanordnung besteht darin, daB mindestens zwolf Dusen 
notwendig sind, um alle Freiheitsgrade gleichzeitig und 
beliebig steuern zu konnen. MuB der Ausfall von zwei beliebigen 
Diisen tolerierbar sein, sind 36 Diisen erf orderlich. 

Demgegeniiber liegt der Vorteil einer geeigneten schiefwink- 
ligen, also nicht orthogonalen und achsenparallelen Dusen- 
anordnung im wesentlichen in folgenden drei Punkten: 

1 . Bei f rei wahlbarer fester Dusenanordnung sind weniger als 

1 2 Diisen erf orderlich . 
2. Fur die Tolerierung des Ausfalls von zwei beliebigen Dusen 
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sind erheblich weniger als 36 Diisen erf orderlich . 
3. Bei treibstof f minimaler Ansteuerung der Diisen ist ein 
geringerer Treibstof fverbrauch moglich. 

Der Nachteil dieser Dusenanordnung besteht darin, dafi die 
treibstof fminimale Dusenans teuerung die Losung eines linearen 
Optimierungsproblems erfordert und der Rechenauf wand fiir die 
Losung dieses Problems mit dem Simplexalgorithmus oft zu hoch 
ist. Aus diesem Grunde sind andere Verfahren entwickelt worden, 
die weniger Rechenzeit benotigen und durch die eine roehr oder 
weniger treibstof fminimierte , aber keine insgesamt treibstof f- 
minimale Dusenans teuerung erreicht werden kann. 

Beiro sogenannten "Linear Programming" wird das Dusenansteue- 
rungsproblem als lineares Optimierungsproblem aufgefafit und mit 
einschlagigen Verfahren der linearen Optimierung, wie zum 
Beispiel den Simplexalgorithmus, geldst. Der schon erwahnte 
Nachteil, der den Einsatz dieses Verfahrens in vielen Fallen 
unmoglich macht, ist der zu hohe Rechenauf wand. 

Beim "Table Look-up" -Verfahren werden alle treibstof fminimalen 
Dusenanordnungen fur Einheitskraf te und -momente bezuglich 
einer Achse des korperfesten Koordinatensystems des Raumfahr- 
zeugs berechnet (Losen des linearen Optimierungsproblems) und 
in einer Tabelle abgelegt. Sind im Einsatz Krafte und Momente 
bezuglich mehrerer Achsen gleichzeitig auszuiiben, werden die 
entsprechenden Diisenansteuerungen bezuglich Kraft Oder Moment 
und fiir jede geforderte Achse und dem geforderten Richtungs- 
Oder Drehsinn aus der Tabelle ausgelesen und mit dem entspre- 
chenden Betragswert der Kraft- bzw. Momentenkomponente in 
dieser Richtung multipliziert . Den Ansteuerwert fur jede Duse 
erhalt man durch die Addition der einzelnen Werte fur die 
jeweilige Duse. Dieses Verfahren ist sehr einfach und sehr 
schnell, fiihrt im allgemeinen aber zu nicht treibstof fminimalen 
Diisenansteuerungen . 

Bei verfeinerten und erweiterten "Table Look-up" -Verfahren 
handelt es sich urn Verfahren, die eine Look-up-Tabelle anlegen, 
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die Eintrage fur zusatzlich eingefuhrte Achsen enthalten. Durch 
diese Verfahren konnen treibstof f sparende Dusenans teuerungen 
ermittelt werden. Zugleich nimmt aber die tibersichtlichkeit des 
Verfahrens ab und wegen der zunehmenden GroBe der Tabellen der 
Rechenauf wand erheblich zu. 

Andere Verfahren beruhen auf der Auswertung und Verbesserung 
des Ergebnisses, das mit dem "Table Look-up" -Verfahren 
ermittelt wird. Das Ziel dieser Auswertung besteht darin, 
besonders treibstof f vers chwendende Dusenans teuerungen zu 
verhindern bzw. geeignet zu korrigieren. Hier soli exemplarisch 
ein Verfahren, das auf einem Puzzy-Ansatz beruht , genannt 
werden (vgl. T. Suzuki , K. Yasuda , S.Yoshikawa, K.Yamada, 
N.Yoshida "An Application of Fuzzy Algorithm to Thruster 
Control Systems of Spacecraft" In: ISTS 94-C-16, 19th Inter- 
national Symposium on Space Technology and Science, Yokohama, 
Japan, May 15-24, 1994, pp. 1-8). 

Eine Treibstof fminimierung der Dusenansteuerung insgesamt hat 
zum Ziel eine solche Ansteuerung der Diisen zu erroitteln, die 

* 

1. die geforderten Krafte und Momente am Raumfahrzeug ausubt, 

2. fur jede Duse einen nichtnegativen Ansteuerwert liefert, 
und 

3. die Summe aller Ansteuerwerte der Diisen minimiert. 

Dies stent - wie erwahnt - ein lineares Optimierungsproblem 
dar. Da das Losen des Problems in der Regel zu zeitauf wendig 
ist, kann die Losung des quadratischen Optimierungsproblems als 
Ersatzproblem durchgefiihrt werden. 

Die Losung des Ersatzproblems fuhrt zur Ermittlung derjenigen 
Dusenansteuerung, die 

1. die geforderten Krafte und Momente am Raumfahrzeug ausiibt, 

2. fur jede Duse einen nichtnegativen Ansteuerwert liefert, 

und 
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3. die Summe aller quadratischen Ansteuerwerte der Diisen 
minimiert . 

Dieses quadratische Optimierungsproblem kann liber die Bereit- 
5 stellung der Moore-Penrose-invertierten Systemraatrix der 

Diisenanordnung am Raumfahrzeug und eines Off setvektors, der mit 
Hilfe der singularen Wertezerlegung berechnet wird (vgl. z-B- 
B.Noble, I,W.Daniel "Applied linear algebra" Prentice-Hall 
International, INC* (1988), pp. 338 - 350), geldst werden. Die 

10 Ldsung des quadratischen Optimierungsproblems kommt dann als 

Produkt der Moor-Penrose-invertierten Systemmatrix rait dem 
kommandierten Kraf te-Moraenten-Vektor und der Addition des 
Of f setvektors zustande. Die Addition des Of f setvektors bewirkt, 
daB im Ldsungsvektor nur nicht negative Ansteuerwerte 

15 auftreten. Die Berechnung des Of f setvektors wird jedoch dann 

numerisch immer aufwendiger, wenn die Zahl der anzus teuernden 
Diisen um zwei oder mehr uber der Dimension des kommandierten 
Krafte-Momenten-Vektors iiegt. Im allgemeinen ist die Losung 
nicht treibstoff -minimal, fallt aber mit der Ldsung des 

20 linearen Optimierungs -problems zusammen, wenn die Zahl der 

anzusteuernden Diisen nur um eins grofier ist als die Dimension 
des kommandierten Krafte-Momenten-Vektors (vergleiche. z.B. 

DE 195 19 732 A1 ) . Weiterhin ist nachteilig, daB langere 
Einschaltzeiten uberproportional "bestraft" werden und daher 

25 seltener in den Ansteuerungen vorkommen. Fiir den Betrieb der 

Diisen und ihrem Wirkungsgrad ist es jedoch giinstiger, sie mit 
langeren Einschaltzeiten zu betreiben. 

Der Nachteil bestehender Verfahren besteht mithin darin, daB 
30 sie entweder nicht das lineare Optimierungsproblem losen und 

damit im allgemeinen nicht die treibs tof f minimale Ansteuerung 
der Diisen bereit s tellen oder aber, wenn bei den Ansteuerver- 
fahren das lineare Optimierungsproblem geldst wird und der 
Simplexalgorithmus ausgefuhrt wird, dieses iiblicherweise zu 
35 rechenaufwendig ist* 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der 
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eingangs angegebenen Art: bereitzustellen, mit dem das lineare 
Optimierungsproblem fur die treibstof fminimale Ansteuerung 
beliebig fester Dusenanordnungen zur gleichzeitigen Reali- 
sierung von Kraften und Moraenten an einem Raumfahrzeug in 
signifikant kurzerer Zeit gelost werden kann, als es bisher 
raoglich ist. 

Diese Aufgfibe wird gemaB der Erfindung durch die im Anspruch 1 
gekennzeichneten Merkmale gelost* 

Das Verfahren nach der Erfindung basiert also vorteilhaft auf 
einem, dem Simplexalgorithmus angelehnten Verfahren, das jedoch 
aufgrund der konsequenten Ausnutzung der Besonderheiten bei der 
Dusenansteuerung bei Raurof ahrzeugen und der Aufgabe der Verwen- 
dung des ublichen Simplextableaus zu einer Rechenzeitverkiirzung 
von ca« einer GroBenordnung fiihrt. 

Vorteilhafte Ausges taltungen und Weiterf uhrungen des Verfahrens 
nach der Erfindung sind in den dem Anspruch 1 nachf olgenden Pa- 
tent anspr lichen gekennzeichnet . 

Das Verfahren nach der Erfindung basiert auf folgenden Grund- 
lagen : 

a) Das bei der Dusenansteuerung zu losende lineare Optimie- 
rungsproblem kann wie folgt auf geschrieben werden: 

Gesucht ist der Ans teuervektor k der die Ans teuerwerte fur 

—opt 

alle Dusen enthalt und folgende Bedingungen erfiillen muB: 



m 



1. Durch JS Qpt miissen die geforderten Krafte und Momente b feo 

am Raumfahrzeug ausgetibt werden. Das heifit, die Glei- 

chung A . k . - b, muB erfiillt werden, wobei A die 

» — sys —opt — kom — sys 

Systemmatrix ist, in der die drei Kraft- und die drei 
Momentenkomponenten fiir jede einzelne Diise in einer Spalte 
dieser Matrix eingetragen sind* Die Systemmatrix hat 6 Zeilen 
wegen der 6 moglichen Freiheitsgrade in der Bewegung des 
Raumf ahrzeugs , namlich 3 Krafte in drei zueinander orthogonalen 
Achsen und entsprechend 3 Momente um diese Achsen und bei 
D Diisen D Spalten. 
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2. k 



Ck 



opt opt 1 * * opt j 

tive Werte enthalten, 

d.h. 



k , ^] darf nur nichtnega- 
opt D * 



* opt j * 0, j - 1(1)D. 
3. die Summe aller Ansteuerwerte muB minimal sein 
( Zielf unlction ) , 



d.h. 



D 
£ k 
j-1 



opt 3 



Min ! 



b) Ein Simplextableau besteht aus einem Spaltenvektor der 

Basisvariablenindizes "X w , einem Zeilenvektor der Nichtbasis- 

variablenindizes % der Koef f izientenmatrix A^ dem 

hd ""tab 

Zeilenvektor der auf die Nichtbasisvariablen urogerechneten 

Zielf unktionskoef f izi^nten J2- tab und dem Spaltenvektor der 

Basisvariablenwerte b^ w : 

•"tab 







— ' hb^+j 






hi 


• 


' * ' *ly 

• 


a lm 

• 






• 

*kl 

• 


• 

a kj 

• 
• 


m 

a km 

• 
• 


• 

b x 


h* 






« 

dm 






P x 


... p 


... p 


-Q 



Im Verfahren nach der Erfindung kommt das Simplextableau nicht 
in dieser vollen Form zum Einsatz. 

c) Fur die Umsetzung eines beliebigen 6-diraensionalen Krafte- 
Momenten-Vektors sind maximal 6 Dusen einer vorhandenen festen 
Dusenanordnung erf orderlich Oder der kommandierte Krafte- 
Momenten-Vektor ist mit der vorhandenen Dusenanordnung nicht 
realisierbar . 

d) Fur jede Dusenauswahl I von 6 aus D Diisen existiert ein 
Basisvariablenindizesvektor "V( I ) . Die Nummern Oder Indizes 

— D 

(1,2,...,D) der in der Auswahl I enthaltenen Diisen werden im 
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Basisvariablenindizesvektor ^ (I) zusammengef afit und die 

b 

Ansteuerwerte dieser Diisen sind die Werte der Basisvariablen 

des Simplextableaus. Die Werte der Basisvariablen sind 

entsprechend der Reihenfolge der Dusenindizes in% b (I) im 

Basisvariablenwertevektor b^ w lV ,b ) abgelegt* Zu jedem 

— tab — d — com 

Basisvariablenindizesvektor *X fo ( X ) gehort eine Basis- 
matrix B b (I). Die Basismatrix wird gebildet, indem man 

diejenigen Spalten der Sys temmat rix A , die zu den 

s y s 

ausgewahlten Diisen gehdren, zu einer neuen Matrix spaltenweise 
zusammens tellt . Die Reihenfolge der Spalten in der Basismatrix 
entspricht der Reihenfolge der Indizes der Diisen im 
dazugehorigen Basisvariablenindizesvektor X b (I). Die 
Nichtbasismatrix B w ( I ) wird aus denjenigen Spalten der 
Systemmatrix zusammengestellt , die zu den nicht in der Auswahl 
I enthaltenen Diisen gehoren* Xhre Indizes sind im Nichtbasis- 
variablenindizesvektor *X . (I) zusammengef afit . Die Reihenfolge 
der Spalten in der Nichtbasismatrix entspricht der Reihenfolge 
der Indizes der Diisen im dazugehorigen Nichtbasisvariablen- 
indizesvektor . Die Ansteuerwerte dieser Diisen sind Null* 

e) 1st fur einen kommandierten Kraf te-Momenten-Vektor ^ Jcom • 

das Optimaltableau * b . opt . ^nb.opt' ^tab.opt' %ab.opt und 
—tab opt 9 efunden > dann sind die Ansteuerwerte der Diisen, die 

durch die Indizes in "X _ ^ adressiert werden, gleich Null und 

— no • opt 

diejenigen Diisen, die den Basisvariablen entsprechen, also 

durch die Indizes in X adressiert werden, enthalten als 

— b • opt 

Ansteuerwert den Wert aus der entsprechenden Zeile des opti- 

malen Basisvariablenwertevektors b^ . . . Der Basisvariablen- 

™"tab ♦ opt 

indizesvektor X. gibt die optimale Diisenauswahl an. 

— b . opt 

f) Fur den Fall der gleichzeitigen Ausiibung von alien drei 
Kraf tkomponenten und alien drei Momentenkomponenten ist der 
Einsatz von maximal 6 der vorhandenen Diisen erforderlich. Es 
gibt genau 



DI 



A(C 6 ) 



6! (D-6)! 



Hoglichkeiten ( Kombinationen ohne Wiederholung) , 6 aus D Diisen 
auszuwahlen. Werden die Indizes in den beiden Vektoren der 
Basisvariablenindizes und Nichtbasisvariablenindizes , *X und 
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\ . , stets aufsteigend sortiert, so gibt es jeweils ge- 
nau A(Cg) dieser Vektoren. Fur eine feste Dusenanordnung gibt 
es (nur) eine bestimmte Anzahl von Kombina tionen von Dusen, die 
die Diisenauswahl fiir eine Losung des linearen Optimierungs- 
problems sein Iconnen. Diese Anzahl ist im aligemeinen deutlich 
Icleiner als die Anzahl AtC^K 

Eine wesentliche Grundlage des Verfahrens besteht darin, daB 
stets mit sortierten Basis- und Nichtbasisvariablenindizes- 
vektoren gearbeitet wird. Im weiteren sind die Indizes Oder 
Nummern der Diisen in den Basis- und Nichtbasisvariablenindizes- 
vektoren z.B. stets der GroBe nach sortiert. Auf diese Weise 
wird es moglich, daB einer Diisenauswahl genau ein Indizes- 
vektor, genau eine Koef f izientenmatr ix A tab ~ ^ a <tab k j ^ 9 genau 
ein Koef f izientenvektpr der auf die Nichtbasisvariablen 
umgerechneten Zielfunktion £ tab und genau eine Basismatrix B fe 
und Nichtbasismatrix B . zugeordnet werden. 

Eine weitere wesentliche Grundlage des Verfahrens besteht 
darin, daB mit einem reduzierten Tableau gearbeitet wird. Das 
reduzierte Tableau besteht lediglich aus der Diisenauswahl- 
identif ikationsnummer I und dem Basisvariablenwertevektor £ tab * 
Fiir das Verfahren und die Arbeit mit dem reduzierten Tableau 
ist neben der Bereits tellung des Feldes der Koef f izienten- 
raatrizen aller Optimaltableaus und der Pivotspaltenindizes- 
matrix auch die Bereits tellung der Dusenauswahlidentif ika- 
tionsnummernmatrix, der Zeilenaustauschindizesmat rix und der 
Basisindizesmatrix erf order lich. 

Das Verfahren besteht aus zwei Bes tandteilen . Der erste, der 
Initialisierungsteil , dient zum einen der Bereits tellung des 
Feldes der Koef f izientenmatrizen aller Optimaltableaus 

£ = 4] = M- wobel Iv 4 = ^ (/) = ■ ** (/) ist - 

Die Dusenauswahlidentif ikationsnummer lauft von 1 bis 

M ( I = 1(1 )M ), wobei M die Anzahl aller Dusenauswahlmoglich- 

keiten ist, die fur einen beliebigen Kraf te-Momenten-Vektor 
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— kom Und die vorlie 9 ende beliebige und unveranderliche 
Diisenanordnung optimal sind. Eine Koeff i2ientenmatrix eines 
Optimaltableaus liegt im Sinne dieses Verfahrens dann vor f wenn 
alle Elemente im zur Koeff izientenmatrix gehorenden Vektor der 
auf die Nichtbas isvariablen umgerechneten 
Zielf unkt ions koeff izient en 

£ tab = ^tab'-tab' nicht negativ sind. Fur die Elemente der 
Koeff izientenmatrix muB daher gelten: 

1 -J/fb M Sp tabj i0fttr ^ 1(116 ' 

2um anderen werden im Initialisierungsteil des Verfahrens 
f olgende Ad refimat risen : die Pivo t spal t enindizesmat rix 
P =» tp T}c ], die Diisenauswahlidentif ikationsnummernma- 
trix U = ^ u ijc^' die Zeilenaustauschindizesmatrix T » Ct Ik ] und 
die Basisindizesmatrix L « ^jjc^ bereitges tellt , damit auf die 
im zweiten Teil des Verfahrens benotigten Daten besonders 
schnell zugegriffen werden kann und bei der Durchfiihrung der 
Basisaustauschschritte die Arbeit mit dem reduzierten Tableau 
moglich wird. Durch die Diisenauswahlidentif ikationsnum- 

mer I ■ 1,2, , M wird jeder Diisenauswahl , die optimal sei-n 

Icann, und damit auch jeder Koeff izieatenmatrix eines Optimal- 
tableaus eine Identif ikationsnummer zugewiesen. 

In der Pivotspaltenindizesmatrix P « CPifc 3 wird fUr 3 ede Zeile 

von jeder in F abgelegten Koeff izientenmatrix der Index der 

dazugehdrigen Pivotspalte s € {0,1,2,3, (D - 6)} abgelegt- 

Fiir eine Koeff izientenmatrix A . und eine gewahlte Pivotzei- 

—tab 

le z ist die Pivotspalte s diejenige, fiir die gilt: 

Minimum P; 

a n — - 



<0 n^T-j^."*/-^^-^ 



D.h. die zu einer Zeile einer Koeff izientenmatrix gehorige 
Pivotspalte ist diejenige Spalte, in der das Mat rixelement der 
Koeff izientenmatrix in dieser Zeile negativ ist und der 
Quotient aus dem Element des Vektors der auf die 

Nichtbasisvariablen umgerechneten Zielf unktionskoef f izienten in 
der Spalte und dem Betrag des entsprechenden Ma trixelements aus 
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der Koef f iz ientenma t r ix am kleinsten ist. 

Mit Hilfe der Pivotspaltenindizesmatrix ist fiir eine Diisen- 

auswahl I, fiir die in £ eindeutig die dazugehorige 

Koef f izientenmatrix 2i tab (I ) abgelegt ist, und gewahlter 

Pivotzeile z die Angabe der entsprechenden Pivotspalte mit 

s(A^ .(I),z) m p aus der Pivot spaltenindizesmatrix moglich, 
""tab Iz 

ohne daC dazu die Koef f izientenmatrix A . w ( I ) Oder der Vektor 

""tab 

der auf die Nichtbasisvariablen umgerechneten Zielfunktions- 
koef f izienten E ta jjt*) bendtigt wird. 

Existiert zu einer bestimmten Zeile keine Pivotspalte, so ist 

das in der Pivotspaltenindizesmatrix durch eine Null (p _ * 0) 

x z 

angegeben. 



In der Ditsenauswahlidentif ikationsnummernmatrix U « Cu ^3 ist 
fiir jede Diisenauswahlidentif ikat ionsnummer I - 1,2, . M und 

jede Zeile z « 1,2,..., 6 jeder, in F abgelegten Koef f izienten- 
matrix die Diisenauswahlidentif ika tionsnummer derjenigen Diisen- 

20 auswahl angegeben, die nach dem Basisaustauschschritt vor- 

liegt, wenn vor dem Basisaustauschschritt die Dusenauswahl I 
vorlag und die Zeile z als Pivotzeile fiir den Basisaustausch- 
schritt gewahlt wurde. O.h., wenn zum Zeitpunkt n die Diisen- 
auswahl I[n3 vorliegt und ein Basisaustauschschritt mit der 

25 Pivotzeile z[n] durchgef iihrt wird, kann die Identif ikations- 

nummer der neu entstehenden Dusenauswahl ohne Verwendung von 

A . (I) mit I[n + 1] - u rn direkt aus der Diisenauswahl- 

— tab I[nj zlnj 

identif ikationsnummernmatrix ausgelesen werden. Dies wird 
dadurch moglich, dafi durch die Kenntnis des aktuellen Tableaus 
30 (bzw. seiner Identif ikationsnummer I) und der Wahl der 

Pivotzeile die Pivotspalte implizit festgelegt ist. Zur 
Ermittlung des Folgetableaus ist also die explizite Kenntnis 
der Pivotspalte s[n] unnotig. 

35 In der Zeilenaustauschindizesmatrix T - ^t Ik 3 ist fiir jede 

Zeile z * 1,2, 6 von jeder der in £ abgelegten Koeffizien- 

tenmatrizen der Zeilenindex des Pivotelements im (sortierten) 
reduzierten Tableau angegeben, das nach dem Basisaustausch- 
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schritt vorliegt, wenn vor dem Basisaus tauschschritt die 
Dusenauswahl I vorlag und die Zeile z als Pivotzeile gewahlt 
wurde. 



In der Basisindizesmatrix L » tij] ist fur jede Dusenauswahl- 
identif ikationsnummer der Vektor der Indizes der in der Dusen- 
auswahl I enthaltenen Dusen, also der zu I gehorige Basis- 
variablenindizesvektor V (I) « 1 _ angegeben. Die Diisenindizes 
in den Vektoren 1^ sind der GrdOe nach sortiert. 

Durch den Initialisierungsteil des Verfahrens wird weiterhin 

eine Start-Dusenauswahlidentif ikat ionsnummer I[0] =» I und 

s tart 

die dazugehorige invertierte Basismatrix 

B w ^4.«^4. 13 < B w(I ^ . >) bereitgestellt. Aus numerischen 
— d • start — b starts 

Griinden 

-1 

werden alle Werte in F und B w . _ # die einen sehr kleinen Be- 

— ^ — d . start 

tragi, (z.B.: £ » 1 0~ ) haben, zu Null gesetzt. 

Iro zweiten Teil des Verfahrens wird zur Realisierung der im . 
Kraf te-Momentenvektor Kommandierten Krafte und Momente die 
optimale Dusenans teuerung , d.h. die Ansteuerwerte der 
einzuschaltenden Dusen ermittelt. Nur dieser Teil des 
Verfahrens rouB den gestellten Echtzeitanf orderungen geniigen. 
Dieser zweite Teil besteht wiederum aus zwei nacheinander 
auszuf uhrenden Teilen. Zuerst wird ein dual zulassiges 
reduziertes Tableau als Starttableau ermittelt. Im zweiten Teil 
wird die erf orderliche Anzahl von Bas isaustauschschritten 
durchgef iihrt . 

Ein dual zulassiges reduziertes Starttableau wird wie folgt 

ermittelt: Fur die Bereitstellung eines dual zulassigen 

reduzierten Tableaus ist lediglich die Transformation des 

aktuellen Krafte-Momente-Vektors b. auf die fur den Start der 

"kom 

Basisaustauschschritte festgelegte Dusenauswahl I[0] = I start 
erf orderlich. Dazu dient die schon berei tgestellte invertierte 
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•1 

Basismatrix B , und es gilt: 

d» stft r t 

-1 

b a B • b 

—tab. start b. start — kom 

Mit der schon bereitgestellten Diisenauswahlidentif ikations- 

nummer ftir den Start ICO] » I und dem neu berechneten Vek- 

start 

tor b. . [03 ah steht ein reduziertes und dual 

—tab — tab. start 

zulassiges Starttableau [I[0] ^ab^ 0 ^ bereit. 

Dann erfolgt die Ausfuhrung der Basisaus tauschschritte (zykli- 
scher Verf ahrensteil ) : 

Der folgend beschriebene Verf ahrensteil muB wiederholt solange 
ausgefiihrt werden, bis die optimale Diisenauswahl mit den 
optimalen Ansteuerwerten gefunden wurde Oder festgestellt wird, 
daB es fur den aktueli vorliegenden Kraf te-Momenten-Vektor und 
die vorliegende Dusenanordnung keine Losung gibt. 

Da von einem dual zulassigen Simplextableau in Form des redu- 
zierten Starttableaus ausgegangen wird und durch die Basisaus- 
tauschschritte immer wieder nur eine Diisenauswahl erzeugt war- 
den kann, deren Simplextableau dual zulassig ist, kann bei der 
Durchfiihrung der Basisaustauschschritte immer wieder auf die im 
Initialisierungsteil des Verfahrens bereitgestellten Daten 
insbesondere im Feld der Koef f izientenmatrizen aller optimalen 
und damit gleichzeitig aller dual zulassigen Simplextableaus 
und in der Pivotspaltenindizesmatrix zuriickgegrif f en werden. 

Erster Schritt: 

Finde einen negativen Wert in b^^Cn]. n in eckigen Klammern 

bezeichnet die Nummer des Iterationsschrittes . Enthalt b^ab'* 1 * 3 

keinen negativen Wert, liegt die Losung des linearen Optimie- 

rungsproblems vor, und es kann der Vektor der Diisenan- 

steuerwerte jL opt aufgebaut werden. Enthalt b^^Cnl negative 

Werte, dann wird eine dieser Zeilen als Pivotzeile z[n] 

gewahlt. Es ist giinstig, die Zeile des betragsgroBten negativen 

Elements in b^ w Cn] als Pivotzeile zu wahlen. 
— tao 
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Zweiter Schritt: 

Mit Kenntnis der Diisenauswahlident if ika t ionsnummer I[n], und 

der im ersten Schritt gewahlten Pivotzeile 2[n] wird aus der 

5 Pivotspaltenindizesmatrix JP der dazugehorige Pivotspaltenindex 

mit sCn] P T r t r t ausgelesen. 1st in P Jcein Pivotspal- 
Xlnjzlnj * — 

tenindex angegeben (p - - P n - 0 ) , mufl eine andere Pivotzeile 

* I[n]z[n] 

ausgewahlt -verden. 1st zu keiner der moglichen Pivotzeiien in 
der Pivotspaltenindizesmatrix eine Pivotspalte indiziert, 
10 existiert keine Losung des linearen Optimierungsproblems . Sind 

Pivotzeile und Pivotspalte ermittelt, wird zum dritten Schritt 
iibergegangen. 

Im dritten und im vierten Schritt wird der eigentliche Basis- 
15 austauschschritt durch Ermittlung des nachsten reduzierten 

Simplextableaus ausgef tihrt : 
Dritter Schritt: 

Mit Hilfe der Koef f izienten in der Pivotspalte s[n] von A , [n] 

wird aus b [n] der neue Basisvariablenwertevektor b . 

— tab "tao.neu 

20 berechnet. Die benotigten Elemente des Pivot spaltenvektors a 

aus ^ tab Cn] Jconnen direkt aus dem Feld der Koef f izienten- 
matrizen aller Optimal tableaus mit a tab fc gj-^Cn] « f fc s[n]I[n] 
ausgelesen werden • 

Die Elemente des neuen Basisvariablenwertevektors werden wie 
25 folgt berechnet: 

Der neue Wert in der Pivotzeile zCn] wird berechnet nach: 

h a btab Z M ^ 

tabneu z[i\\ 7 — — , 

Jz[n\s[n\I[n) 

30 wobei f r n r ^ ^ r . das entsprechend indizierte Pivotelement 

z[nj sCn] I[n] * 

au s . Cn] ist. 

"tab 

Die neuen Werte in den anderen Zeilen werden berechnet nach: 
3 5 b tabncu k = b tab k M ~ b tabjieu z[n] ' fks[n]I{n] 

wobei: (Ac — 1(1)6) aber [k * z[n]) 

Da der neuberechnete Basisvariablenwertevektor _ M zum 

~tab.neu 
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fiktiv angenommenen neuen sortierten Basisvariablen- 
indizesvektor ^. b Cn + 13 passen muG, ist im allgemeinen der 
neuberechnete Basisvariablenwertevektor entsprechend 
umzuordnen* Dazu sind im neuberechneten Basisvariablen- 
wertevektor b,_ . lediglich die Werte zwischen der aktuellen 
Pivotzeile z[n3 und der mit Hilfe von I[n] und zCn] aus der 

Zeilenaustauschindizesmatrix T - [t T 3 auszulesenden Austausch- 

— I z 

zeile z^_ . . Cn] « t ^ r , r zyklisch zu tauschen. Die 
tau£ch I[n] z[n] * 

Austauschzeile z^ ^[n] ist die Zeile, wohin der Wert des 

tausch 

neuberechneten Basisvariablenwertevektors plaziert wird, der in 

b x . in der Pivotzeile zCn3 steht. Nach diesem Umordnen des 

—tab . neu 

neuberechneten Basisvariablenwertevektors liegt der fiir den 
f olgenden Basis aus tauschschritt benotigte Basisvariablen- 
wertevektor b. w [n + 1 ] vor. 

Vierter Schritt: 

Mit Kenntnis von Itn] und der Pivotzeile z£n] wird mit 

I[n + 1] - u Tr n r . die Diisenauswahlidentif ikationsnummer 
XLnj zlnj 

ermittelt, die zum neuberechneten und nach aufs teigenden 
Basisvariablenindizes sortierten Basisvariablenwertevek- 
tor £l tab E n + 13 gehort, Der Basisaustauschschritt ist damit . 
vollzogen, das neue reduzierte Tableau [I[n +13 ^ tab t n + 1 ^ 
liegt vor, und es kann die Losung des linearen Optimierungs- 
problems mit dem ersten Schritt des zyklischen Verf ahrensteils 
fortgesetzt werden. 

Bricht das Verfahren mit einem Basisvariablenwertevektor, der 
keine negativen Werte mehr enthalt, ab, so liegt der optimale 
6-dimensionale Basisvariablenwertevektor li tat)#0 p t " -tab^ n end^ 
vor, und es kann mit Hilfe von I C n end ^ aus der Basisindizes- 
matrix L der zum aktuellen Basisvariablenwertevektor gehorende 
Vektor der Basisvariablenindizes mit *^b[n d l = 1 I [ n ] 

end 

ausgelesen werden. Mit Kenntnis der beiden Vektoren *i tab E n enci ^ 

und "X w [n ^3 wird der optimale D-dimensionale Ans teuervektor 
— d end 

fiir die Diisen j£ opt aufgebaut. 

Fiir die Losung des linearen Optimierungsproblems fiir 
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unterschiedliche Krafte-Momenten-Vektoren sind im allge- 

meinen auch unterschiedlich viele Basisaustauschschrit- 

te n . € {0, 1 , . .., n } erf order lich . Bei Simulationsrech- 
ena max 

nungen mit einer Dusenanordnung von 12 Diisen vrurden in sehr 

seltenen Fallen bis zu n = 10 Basisaus tauschschritte durch- 

max 

gefuhrt, obwohl im Durchschnitt nur ca. 3,5 Basisaus tausch- 
schritte notig waren. Bei geeigneter Wahl des Starttableaus 

CI ^ b . . . A ^ ] sind jedoch nur in ca. 0,5 % aller Falle 

start —taD, start 

mehr als 5 Austauschschritte notwendig. Im allgemeinen bleibt 
es ohne roerUiche Konsequenzen fiir die Anwendung, wenn die 
Anzahl der maximal auszuf iihrenden Austauschschritte auf einen 
Wert begrenzt wird, bei dera nahezu alle Falle vollstandig 
bearbeitet verden und in den Fallen, in denen doch noch 
negative Werte im Basisvariablenwertevektor stehen, diese zu 
Null gesetzt werden. Auf diese Weise Itann der Worst-Case- 
Rechenaufwand weiter gesenkt werden. Anwendungen dieser 
Verf ahrensmoglichkeit bei der Simulation in geschlossenen 
Regelkreisen zeigten keine bemerkbare Verschlechterung der 
Regelgiite oder eine Erhohung des Treibstof f verbrauchs . 

Enthalt der Basisvariablenwertevektor nach der Ausfuhrung der 
maximal erlaubten Anzahl von Basisautauschschritten noch 
negative Werte, ist es prinzipiell auch moglich, daO keine 
Losung des linearen Optimierungsproblems existiert. Dieser Fall 
kann dann auftreten, wenn eine Dusenanordnung vorliegt, die 
nicht beliebige Krafte und Momente aus dero Bereich aller 
zugelassenenen Kraf te-Momenten-Vektoren erzeugen kann. Damit in 
diesen Fallen keine zu groBen Fehlkrafte und Fehlmomente 
erzeugt werden, kann in Abhangigkeit vom Betrag der negtiven 
Elemente im Basisvariablenwertevektor entschieden werden, ob 
die negativen Elemente zu Null gesetzt werden und der Ansteuer- 
vektor fur die Diisen aufgebaut wird oder generell keine An- 
steuerung der Diisen erfolgt. 

Kommen bei der Abarbeitung des Simplexalgorithmus entartete 
Optimaltableaus vor, dann ist es theoretisch moglich, daB der 
Algorithrous in einen Endloszyklus lauft und nicht konvergiert. 
In Simulationsrechnungen mit unterschiedlichen Diisenanordnungen 
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ist dieser Fall allerdings noch nie eingetreten und wird 
deshalb auch nicht explizit beriicksichtigt . Fur den dualen 
Simplexalgorithmus kann man durch die gewahlte Dusenanordnung 
sicherstellen, daO keine dual entarteten Optimaltableaus 
existieren. Fiir den primalen Simplexalgorithmus, also das 
Auftreten primal entarteter Optimaltableaus, ist das so nicht 
raoglich. 

Soil bei der Ermittlung der treibstof fminimalen Triebwerks- 
ansteuerung der unterschiedliche spezifische Treibstoff- 
verbrauch c ^ , j»1(1)D der einzelnen Triebwerke beriicksichtigt 
werden, liegt die Zielfunktion 

vor • J 

Die Berucksichtigung der Gewichtskoef f izienten ist in dem 
Verfahren £ur den on-line auszuf iihrenden Verf ahrensteil ohne 
Mehrauf wand moglich. 

Fallt eine Diise aus , andert sich die Dusenanordnung im Sinne 
dieses Verfahrens. Das Verfahren auf die neue Dusenanordnung 
anzuwenden erfordert, daB alle Daten, die im Initiali- 
sierungsteil des Verfahrens bereitgestellt wurden, fur die neue 
Dusenanordnung bereitgestellt werden. 

Da sich die neue Dusenanordnung nur durch den Wegfall einer 
Diise von der alten Dusenanordnung unterscheidet, gibt es 
Moglichkeiten, die neuen Daten durch eine Aktualisierung der 
alten Daten zu ermitteln. Durch die Bereitstellung eines 
erweiterten Feldes von Koef f izientenmatrizen, in dem auch 
diejenigen Koef f izientenmatrizen abgelegt sind, die durch den 
Ausfall von Dusen Koef f izientenmatrizen von Optimaltableaus 
werden, und entsprechend erweiterter AdreBmatrizen, kann der 
Aktualisierungsaufwand stark reduziert werden und fur bestimmte 
Matrizen ganz entfallen. 



SchlieBlich besteht die Mdglichkeit, die Daten fiir alle 
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Diisenanordnungen, die durch den Ausfall einer Diise zu erwarten 
sind, vorab bereitzustellen , und bei Ausfall einer bestimmten 
Diise dann auf die entsprechenden neuen Daten umzuschalten . Bis 
zum Ausfall einer nachsten Diise miissen fur die dann zu 
erwartenden Diisenanordnungen die Daten wiederum bereitgestellt 
werden, urn dann durch Umschalten auf den nachsten Dusenausfall 
reagieren zu konnen. 

Auch bei einer Verschiebung des Schwerpunktes des Raumfahr- 
zeuges kommt es zu einer Veranderung der Dusenanordnung im 
Sinne dieses Verfahrens. Im allgemeinen ftthrt das zu einer 
Veranderung der Systemmatrix. Daher miissen alle von ihr abge- 
leiteten Daten neu bereitgestellt werden. 

Das zuvor beschriebene Verfahren kann zur weiteren Einsparung 
von Rechenaufwand und Speicherplatzbedarf noch modifiziert 
werden : 

- In voranstehender Beschreibung wurde stets von der Steuerung 
aller 6 Freiheit sgrade des Rauraf ahrzeugs ausgegangen. Das 
Verfahren ist auch auf die Steuerung von 1 bis 5 Freiheits- „ 
graden anwendbar. Es ist lediglich die Streichung bzw. die 
Nichtberucksichtigung der den nicht zu steuernden Freiheits- 
graden entsprechenden Zeilen in der Systemmatrix und der von 
ihr abgeleiteten Vektoren und Matrizen erf orderlich . 
Desgleichen ist die Steuerung einer sich standig andernden 
Teilmenge aller 6 Freiheitsgrade moglich. Liegt eine schiefe 
Dusenanordnung vor, bleiben jedoch die nicht gesteuerten 
Freiheitsgrade im allgemeinen nicht unbeeinf luBt . 

- Miissen am Raumfahrzeug nicht beliebige Krafte und Momente 
ausgeiibt werden, ist der Einsatz einer reduzierten Dusen- 
anordnung moglich. Das vorliegende Verfahren ist auch fur 
derartige Anordnungen anwendbar. Es sollte nur gewahrleistet 
sein, daB die kommandierten Krafte und Momente des Kommandos 
durch die vorliegende Dusenanordnung realisierbar sind. Ist da 
Kommando nicht realisierbar, gibt das Verfahren an, dafl keine 
Losung des linearen Opt imierungsproblems existiert. 
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- In Abhangigkeit von der vorliegenden Diisenanordnung kommt es 
vor, daO von den im Feld der Koef f izientenmatrizen aller 
Optima It ableaus abgelegten Koef f izientenmatrizen bestiramte 
Spalten nie als Pivotspalte ausgewahlt werden. Nach Berechnung 
der Pivotspaltenindizesmatrix P konnen die nicht benotigten 
Spalten aller Koef f izientenmatrizen einfach ermittelt werden. 
Es ist moglich, ein Feld aller relvanten Pivotspaltenvektoren V 
zu erstellen, in dem nur die Spalten der Koef f izientenmatrizen 
aller Optimaltableaus abgelegt sind, die tatsachlich als Pivot- 
spaltenvelctoren benotigt werden. Dadurch kann eine signif iJcante 
Speicherplatzersparnis erreicht werden. Die Adressierung der 
Matrixspalten im Feld aller relevanten Pivotspaltenvelctoren 
erfordert keinen erhdhten Rechenauf wand. Jede Spalte hat im 
Feld aller relevan-ten Pivotspaltenvelctoren ihren Spaltenindex. 
Wird in einer veranderten Pivotspaltenindizesmatrix, der 
Pivotspalteniden-tif ikationsnummernmatrix 0, = E°ik^ nicht der 
Spaltenindex eines Pivotspaltenvektors in einer 
Koef f izientenmatrix angegeben, sondern gleich die 
Pivotspaltenidentif ikat ionsnummer als Spaltenindex r im Feld 
aller relevanten Pivotspaltenvektoren V, dann kann im zweiten 
Schritt des zyklischen Verf ahrensteils mit r = o Jk gleich der 
Spaltenindex im Feld aller relevanten Pivot-spaltenvektoren V * 
Cv ] angegeben werden. Die Ermittlung der Elemente des neuen 
Basisvariablenwertevektors im dritten Schritt des zyklischen 
Verf ahrensteils andert sich dann wie folgt: 

, _ btab f&ll ^ 

tabsi*uz[n] V^j^j » 

wobei v r r , das entsprechend indizierte Pivotelement aus 

z[n] rtnj ^ 

A^ w Cn] ist. Die neuen Werte in den anderen Zeilen werden 
"tab 

berechnet nach: 

b tabjtcu k ~ b tab * I/O - b tabjtai z\n\ * V k f{n\ 

wobei : (k « 1(1)6) aber [k * z[n]) 
Die weiteren Verf ahrensbestandteile bleiben unverandert. 

- Liegt ein Feld aller relvanten Pivotspaltenvektoren V vor, 
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Icann insbesondere bei Diisenanordnungen, die bestimmte 
Symroetrien aufweisen, festgestellt werden, dan relevante 
Pivots pal tenvektoren unterschiedlicher Optima Itableaus 
numerisch identisch sind. Da das mehrmalige Ablegen gleicher 
Pivotspaltenvektoren nicht notwendig ist, kann erneut weiterer 
Speicherplatz eingespart werden. Beim sukzessiven Aufbau des 
Feldes aller relevanten Pivotspal tenvektoren im Initiali- 
sierungsteil des Verfahrens kann uberpriift werden, ob ein 
numerisch identischer Pivotspaltenvektor zura gerade ermittelten 
Pivotspaltenvektor schon in V abgeiegt wurde, 1st dies der 
Fall, kann das erneute Ablegen dieses relevanten 
Pivotspaltenvektors entfallen und in der zu V gehorigen 
Pivotspaltenidentif ikationsnummernmatrix £ wird die 
Pivotspaltenindentif ikationsnummer ( Spal tenindex ) des schon im 
Feld aller relevanten Pivotspaltenvektoren abgelegten 
Pivotspaltenvektors eingetragen. 

- Besteht das Erfordernis, den Speicherbedarf des Verfahrens 
noch weiter zu reduzieren, kann eine Verf ahrensvariante 
eingesetzt werden, bei der es nicht erforderlich ist, alle 
relevanten Pivotspaltenvektoren, sondern nur noch Reprasen- > 
tantenvektoren im Rechnerspeicher abzulegen. ReprMsentanten- 
vektoren sind relevante Pivotspaltenvektoren, aus denen mit nur 
minimalem Rechenauf wand weitere relevante Pivotspaltenvektoren 
berechnet werden konnen. 

T 

Jeder Pivotspaltenvektor a = Ca i a 2 *"* a i a n^ 9 eh5rt zu 

einer bestimmten Basis, z.B. 

W ' CX k. 1*1.2 •••^p.n ]T Und 2U einer 

Nichtbasisvariablen, z.B.X w = - Wenn a ^ verschieden von 

no r i rj> 

Null ist, laBt sich aus ad^^) - C a 1 a 2 ** # a i ### a n 3 
derjenige Pivotspaltenvektor 

T 

a = a(I *K ) = [a a ~ ••• a ••• a ] mit 

^neu 2. KX 2* s La neu.1 neu.2 neu.x neu.n 



a neui = Zf: und a^jjjc * i) - - a k • a^j 

berechnen, der zur Basis 
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Nichtbasisvariablen "X- . «X gehdrt. So kann aus a(I-, "X ) mit 

no s — • i r 

der Angabe des Austauschindexes i der neue Pivotspaltenvektor 
all,, "X ) berechnet werden. 

Auf diese Weise konnen aus einero z.B. sechsdimensionalen 
Pivotspaltenvektor bis zu sechs weitere Pivotspaltenvektoren 
berechnet werden* Ftir die Anwendung im Verfahren zur treib- 
stoff minimalen Diisenans teuerung kann man sich jetzt auf das 
Speichern derjenigen relevanten Pivotspaltenvektoren als 
Reprasentantenvektoren beschranken, die erforderlich sind, um 
alle relevanten Pivotspaltenvektoren wie angegeben zu 
berechnen. 

Fiir die Ermittlung einer mdglichst kleinen Anzahl von 
Reprasentantenvektoren lassen sich sehr giinstig Genetische 
Algorithmen als Optiomierungsverf ahren einsetzen. 

Fiir eine geeignete Adressierung der Reprasentantenvektoren in 
der Matrix aller Repra sent anten sind statt der 
Pivotspaltenindizes - bzw. 

Pivotspaltenindentif ikationsnummernmatrix einmal eine 
Reprasentantenidentif ikationsnummernmatrix R und weiterhin eine 
Aus tauschindizesmatrix S, erforderlich. In der 
Reprasentantenidentif ikationsnummernmatrix R wird fiir jede 
Diisenauswahlidentif ikationsnummer und jede Pivotzeile die 
Identif ikationsnummer des Reprasentantenvektor s angegeben, der 
fiir die Berechnung des tatsachlich benotigten 
Pivotspaltenvektors verwendet werden kann* In der 
Austauschindizesmatrix S wird fiir jede 

Diisenauswahlidentif ikationsnummer und jede Pivotzeile der 
Zeilenindex des Aus tauschelements angegeben, um aus dem 
Reprasentantenvektor den tatsachlich benotigten Pivotspalten- 
vektor berechnen zu konnen* Steht an einer Stelle in der 
Austauschindexmatrix eine Null, bedeutet dies, daB der 
adressierte Reprasentantenvektor unmittelbar der erf orderliche 
Pivotspaltenvektor ist und daher die Umrechnung in den aktuell 
benotigten Pivotspaltenvektor entfallt . 
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Da stets rait nach Basisvariablenindizes sortierten 
Pivotspaltenvelttoren gearbeitet wird, muB der zunachst aus dera 
Reprasentantenvektor berechnete Pivotspaltenvektor noch 
umsortiert werden (z.B. nach auf steigenden 

Basisvariablenindizes), urn benutzt werden zu konnen. Wird die 
Sortierung der GroBe der Basisvariablen vorgenoramen , gibt es 
bei einera sechselementigen Vektor genau 30 Moglichkeiten der 
Umsortierung, um die erf or-derliche Ordnung wieder 
herzustellen. Bei einer erforderlichen Umsortierung wird stets 
ein Vektorelement an einer Position im Vektor entnomroen und auf 
eine andere Position geschrieben. Die Vektoreleraente zwischen 
den beiden Positionen werden entspre-chend verschoben. 

In einer Sortieridentif ikationsnummernmatrix G « [g z ^3 wird 
fur jeden Repra sentanten und jede Austauschzeile die Sortier- 
identif ikationsnummer angegeben, nach der der Vektor umgeordnet 
werden muB, urn den benotigten Pivotspaltenvektor zu erhalten. 

Aus einera Reprasentantenvektor konnen durch Basisaustausch 
und/oder oben beschriebene Umordnung bis zu 7 x 31 =217 
Vektoren gebildet werden* Werden alle diese Moglichkeiten 
berucksichtigt , kann es aufgrund von Symmetrien in der 
Triebwerksanordnung vorkommen, daB weniger Reprasentanten- 
vektoren erforderlich sind, als wurden nur die sieben aus einem 
Reprasentantenvektor gebildeten relevanten Pivotspaltenvektoren 
in Betracht gezogen. 

Eine weitere Verf ahrensvariante ergibt sich daraus, daB auf den 
Basisaustausch am Reprasentantenvektor verzichtet wird und nur 
die 30 Umordnungsmoglichkeiten des Reprasentantenvektors 
beriicksichtigt werden. Auf diese Weise reprasentiert ein 
Reprasentantenvektor 31 Vektoren. Aufgrund von Symmetrien in 
der Triebwerksanordnung konnen auch auf diese Weise mehrere 
relevante Pivotspaltenvektoren generiert werden. Im Vergleich 
zur letzten Verf ahrensvariante werden im allgemeinen jedoch 
raehr Reprasentantenvektoren benotigt, jedoch entfallt 
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gleichzeitig der Rechenauf wand fur den Basisaus tausch am 
Reprasentantenvektor r und es ist die Aus tauschindizesmatrix £ 
nicht erf orderlich. 



RNsnnnin- <wo 9BA9058Ai i > 
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Pa t en t an a prtiche 

1 . Verf ahren zur insgesamt treibstof fminimalen rechnerge- 

stiitzten Ansteuerung einer beliebig festgelegten Anzahl und 

beliebig festgelegten Anordnung von Dusen an einem 

Raumf ahrzeug , 

dadurch gekennzeichnet , 

- dafi ein treibstof fminiraaler Ansteuervektor ( k . ) mit 

opt 

ausschlieBlich nicht negativen Ansteuerwerten fur alle 
vorgegebenen Dusen mit Hilfe des dualen Simplexalgorithmus 
gebildet wird, 

- dafi dazu jeder Dusenauswahl aus der Gesamtzahl von Dusen am 

Raumf ahrzeug, die die optimale Dusenauswahl fiir die treib- 
stof fminimale Realisierung eines beliebigen Kommando- 
vektors (b w ) darstellt, eine Diisenauswahliden- 
tif ikationsnummer (X) zugeordnet wird, 

- daB fiir die Ermittlung der optimalen Ansteuerung zunachst 

der Kommandovektor (b, ) der die auszuubende Kraft 

— kom 

und/oder Momentenkomponenten enthalt, durch Multiplikation 

mit einer Transf ormationsmatrix (<B . ) in einen 

start 

Startvektor (b=b \ transf ormiert wird, 

— — start 

- daB danach der duale Simplexalgorithmus durchgefuhrt wird, 

bis der Vektor der Ansteuerwerte (b) nur nichtnegative 
Ansteuerwerte enthalt Oder aber die Nichtexis tenz einer 
Losung angezeigt wird, 

- und daB dabei der duale Simplexalgorithmus realisiert wird, 

indem vorab alle Daten, die fiir die Abarbeitung des dualen 
Simplexalgorithmus notwendig sind und sich off-line (vor- 
ab) berechnen lassen, ermittelt werden, in Datenfeldern 
im Rechnerspeicher abgelegt werden und beim Einsatz des 
dualen Simplexalgorithmus wahrend des on-line-Betriebes 
durch geeignete Adressierungen aus dem Rechnerspeicher 
abgerufen und in der Rechnung verwendet werden, wobei nur 
mit einem reduzierten Simplextableau gearbeitet wird, in- 
dem lediglich die Dusenauswahlidentif ikationsnummer und 
der dazugehorige Vektor der Ansteuerwerte (b) gefuhrt 
werden. 
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2* Verfahren nach Anspruch 1 , 
dadurch gekennzeichnet, 

da0 die Transf ormat ionsmat rix (<B . die invertierte 

start 

Basisroatrix derjenigen DUsen darstellt, die zur Startdiisen- 

auswahl gehdren, diese Diisenauswahl fur einen beliebigen 

Kommandovektor optimal ist, diese Diisenauswahl eine Start- 

Diisenauswahlidentif ikationsnummer < I start ) zugeordnet ist und 

der Startvektor ( b*b - > die Ans teuerwerte dieser Dusen 

— ""Start 

enthalt. 



3. Verfahren nach einem der Anspruche 1 Oder 2, 
dadurch gekennzeichnet , 

- dafi jede Diise am Raumfahrzeug eine Identif ikationsnummer 

erhalt, und die Werte im Vektor der Ans teuerwerte fiir eine 
bestimmte Diisenauswahl stets in einer festgelegten Ordnung 
anhand der Identif ikationsnummern der Diisen sortiert 
werden, 

- daB bei der Abarbeitung des dualen Simplexalgorithmus nach 

der Ermittlung einer Pivotzeilennumraer und mit Kenntnis 
der Diisenauswahlidentif ikationsnummer der aktuell vor- 
liegenden Diisenauswahl die entsprechende Pivotspalten- * 
nummer aus einer vorab bereitgestellten Pivotspaltenin- 
dizesmatrix ausgelesen wird, 

- daft das Nichtvorhandensein einer Pivotspaltennummer zu der 

vorliegenden Diisenauswahlidentif ikationsnummer und zu 
der gewahlten Pivotzeile durch einen besonderen Eintrag 
in der Pivotspaltenindizesmatrix eindeutig angezeigt wird, 

- und daB die Nichtexistenz einer Losung dadurch angezeigt 

wird, daB zu keiner moglichen Pivotzeile in der Pivotzei- 
lenindizesmatrix eine Pivotspaltennummer angegeben ist. 

4. Verfahren nach einem der Anspruche 1 bis 3, 
dadurch gekennzeichnet , 

- daB der Pivot spaltenvektor der Koef f izientenmatrix des 

aktuellen Simplextableaus , dessen Werte fur die Berechnung 
neuer Ansteuerwerte fiir die Diisen einer neuen Diisenauswahl 
bendtigt werden, mit Kenntnis der Diisenauswahl vor dem 
Basisaustauschschritt und der gewahlten Pivotspaltennummer 
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aus einem vorab bereitges tellten Feld aller moglichen 
optimalen Koeffizientenraatrizen ausgelesen wird, 

- daB der Vektor der Ansteuerwerte fur die Dusen der Dusen- 

auswahl, die nach dem Basisaustauschschritt vorliegt, 
berechnet wird und die Werte in diesem Vektor entsprechend 
der Zuordnung zu den Dusen und entsprechend der festgeleg- 
ten Ordnung der Identif ikationsnummern fiir die Diisen 
sortiert werden. 

5. Verfahren nach einem der Anspriiche 1 bis 4, 
dadurch gekennzeichnet , 

- daB mit Kenntnis der Diisenauswahlidentif ikationsnummer vor 

dem Basisaustauschschritt und der gewShlten Pivotzeilen- 
nummer die Diisenauswahlidentif ikationsnummer derjenigen 
Diisenauswahl, die nach dem Basisaustauschschritt 
vorliegt, aus einer Diisenauswahlidentif ikationsnummern- 
matrix ausgelesen wird 

- und daB nach Vorliegen eines Vektors der Ansteuerwerte, der 

nur noch nicht negative Werte enthSlt, mit Kenntnis der 
dazugehSrigen Diisenauswahlidentif ikationsnummer der Vektor 
der Dusenindizes aus einer Basisvariablenindizesmatrix • 
ausgelesen wird und die mit ihrer Dusenidentif ikations- 
nummer indizierten Diisen den entsprechenden Wert aus dem 
Vektor der Ansteuerwerte (b) erhalten sowie die ubrigen 
Dusen den Ansteuerwert Null zugewiesen bekommen und damit 

der Ansteuervektor ( k . ) bestimmt ist. 

~~opt 

6. Verfahren nach einem der Anspriiche 1 bis 5, 
dadurch gekennzeichnet , 

daB nur diejenigen Spalten von Koef f izientenmatrizen dual 
zulassiger Simplextableaus in einer Matrix aller relevanten 
Pivotspaltenvektoren abgelegt werden, die bei der Abarbeitung 
des Verfahrens tatsMchlich als Pivotspaltenvektoren ausgewahlt 
werden kdnnen, 

und daB identische relevante Pivotspaltenvektoren nur einmal in 
dieser Matrix abgelegt werden* 
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7. Verfahren nach Anspruch 6, 
dadurch gekennzeichnet , 

- daB fiir die Adressierung der relevanten Pivotspalten- 

vektoren jedem Pivotspaltenvektor in der Matrix aller 
relevanten Pivotspal tenvektoren eine Pivotspalteniden- 
tif ikationsnummer zugewiesen wird und diese in einer 
Pivotspaltenidentif ikationsnummernmatrix derart eingetra- 
gen werden, daB bei Kenntnis der aktuellen Dusenauswahl- 
identif ikationsnummer und der aktuell gewShlten Pivotzei- 
lennummer die Pivotspaltenidentif ikationsnummer des aktu- 
ell bendtigten Pivotspaltenvektors auslesbar ist, 

- dafi das Nichtvorhandensein einer Pivotspalte zu einer 

gewahlten Pivotzeile durch einen besonderen Eintrag in der 
Pivotspaltenidentif ikationsnummernmatrix eindeutig 
angezeigt wird, 

- und daB die Nichtexistenz einer Losung dadurch bestimmt 

wird, daB zu keiner moglichen Pivotzeile in der Pivot- 
spaltenidentif ikationsnummernmatrix eine Pivotspalten- 
identif ikationsnummer angegeben ist . 

8. Verfahren nach einem der Anspruche 3 bis 7 , 
dadurch gekennzeichnet , 

daB statt der relevanten Pivotspaltenvektoren nur bestimmte 
Pivotspaltenvektoren als Reprasentantenvektoren im 
Rechnerspeicher abgelegt werden, wobei die 
Reprasentantenvektoren eine moglichst kleine Anzahl von 
Pivotspaltenvektoren darstellen, aus denen die tatsachlich 
benotigten relevanten Pivotspaltenvektoren ermittelt werden, 
und daB jeder Reprasentantenvektor eine 

Reprasentantenidentif ikationsnummer zugewiesen bekommt und zur 
Adressierung der Reprasentantenvektoren eine 

Reprasentantenidentif ikationsnummernmatrix verwendet wird, in 
der fiir jede Dusenauswahlidentif ikationsnummer und jede 
Pivotzeilennummer die Reprasentantenidentif ikationsnummer 
desjenigen Reprasentantenvektors eingetragen ist, der fiir die 
Berechnung des erf orderlichen Pivotspaltenvektors verwendet 
wird. 
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9. Verfahren nach Anspruch 8, 
dadurch gekennzeichnet, 

daB in einer Austauschmatrix fiir jede Diisenauswahl- 
identif ikationsnummer und jede Pivotzeilennummer der 
Zeilenindex der auszutauschenden Basisvariablen als 
Austauschelement angegeben ist, um aus dem ReprasentantenveUor 
den benotigten Pivotspaltenvektor zu berechnen. 

10. Verfahren nach einem der Anspruche 1 bis 9, 
dadurch gekennzeichnet # 

daB fttr die Moglichkeit der Veranderung der Diisenanordnung 
durch zu erwartende Dtisenausf alle die entsprechenden Daten 
vorab berechnet und gespeichert werden und nach Feststellen 
eines Dtisenausf alls automatisch durch den Rechner auf Verwen- 
dung der entsprechenden neuen Da-ten umgeschaltet wird. 

11. Verfahren nach Anspruch 10 f 
dadurch gekennxeichnet , 

daB die nach einem Dtisenausf all erf orderlichen neuen Daten 
durch eine entsprechende KorreKtur und Aktualisierung der bis 
zum Diisenausfall verwendeten Daten bereitgestell t werden. 

12. Verfahren nach einem der Anspruche 1 bis 11, 
dadurch gekennzeichnet , 

daB das Verfahren unabhangig von der Dimension des Kommando- 
vektors (b w ) angewendet wird und dabei fur die Berechnung der 
vorab bereitgestellten Daten die nichtauszuUbenden Komponenten 
aus der vollstandigen Systemmatrix der Diisenanordnung 
eliminiert werden. 
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